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ABSTRACT 

We provide a unified description of thermal equilibria of black hole accretion disks, includ- 
ing the newly-discovered advection-dominated solutions. We classify the solutions on the 
basis of optical depth and importance of advection cooling. We demonstrate that only four 
physically distinct topological types of equilibria exist. Two of the types correspond to 
optically thin and optically thick equilibria, while the other two types are distinguished by 
whether advection is negligible or dominant. A stable Shakura-Sunyaev disk exists only 
for accretion rates M below a certain maximum. However, there is a critical viscosity 
parameter a C rit, which is a function of radius, such that for a > a C rit advection-dominated 
solutions exist for all M. Even when a < a C rit, the advection-dominated solutions are 
available for a wide range of M except for a gap around the Eddington rate. We therefore 
suggest that advection-dominated flows may be more common than standard thin disks 
in black hole systems. For certain ranges of radii and M, no stable steady state solution 
is possible. In these cases, we suggest that limit cycle behavior may occur, leading to 
variability. 

Subject headings: accretion, accretion disks — instabilities 
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1. INTRODUCTION 



Accretion disks are common in a variety of astrophysical systems, including young stars, 
cataclysmic variables, X-ray binaries, and active galactic nuclei. A simple model of an 
accretion disk was developed by Shakura & Sunyaev (1973) based on the assumption that 
the energy released through viscous dissipation is radiated locally and efficiently, so that 
the disk vertical thickness is much smaller than the radius. This "thin disk" model has 
been the mainstay of accretion studies for many years. However, the model is known to 
suffer from thermal and viscous instabilities under certain circumstances (Pringle et al. 
1973, Lightman & Eardley 1974). Shapiro, Lightman & Eardley (1976) suggested that 
these instabilities sometimes result in a two-temperature disk. 

Starting with the initial work of Abramowicz et al. (1988) and the more recent 
papers of Narayan & Popham (1993), Narayan & Yi (1994, 1995a) and Abramowicz et 
al. (1995), it has been demonstrated that an unstable accretion disk may switch to a 
completely different mode of accretion where the viscously released energy is advected 
with the gas rather than being radiated. The resulting "advection-dominated" flow is 
stable (Abramowicz et al. 1995, Narayan & Yi 1995b) both to the thermal and viscous 
instabilities. Advection-dominated flows are much hotter than cooling-dominated thin 
disks and are therefore thicker in the vertical direction. For this reason they are also 
described as "thick disks" or "slim disks," depending on the relative vertical thickness. 

Since we now have at least two distinct stable accretion solutions, namely cooling- 
dominated flows and advection-dominated flows, it is important to map out which kind 
of solution occurs under different conditions. This is the purpose of the present Letter. 
We demonstrate the topological relationships among the various solution branches in the 
case of accretion onto a black hole, and unify previous discussions by including the new 
advection-dominated solutions. 



We use the slim disk formalism described in Abramowicz et al. (1995). This model makes 
the following assumptions: (i) The vertical half-thickness of the disk H is smaller than 



calculated by using a pseudo-Newtonian potential (Paczyhski & Wiita 1980). (ii) The 
kinematic viscosity coefficient is written as v = (2/3)ac s H . (iii) The disk is assumed to 
rotate at O = Qk- 



2. BASIC EQUATIONS 






We use the following three steady state equations: 

(angular momentum conservation) (2-1) 

where S = 2Hp is the surface density, /* = 1 - 90*/(0(i?/i? G ) 2 ) with R G = 2GM/c 2 , 
g = — (2/3)(<fln fi/c/ln R) } and fi* is equal to fi#- at i? = 3i?G! 

M = — 27ri?Eu r , (mass conservation) (2-2) 

where u r is the radial velocity; and 

= Q_ + Qadv, (energy conservation) (2-3) 

where the left hand side is the viscous heating per unit area and the right hand side is the 
sum of radiative cooling and advective cooling. In conventional cooling-dominated thin 
disks Qadv is negligible, but this term dominates in advection-dominated flows. We take 

Q + = — &f*g, Q^ v = -Vv r % = (2.4) 
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where £ = -[(4 - 3/3)/(r 3 - 1 )](d In T/d In R) + (4 - 3/3)(</ln E/c/lni?), T is the disk mid- 
plane temperature, f3 = p ga s/.P is the ratio of gas to total pressure, T3 = 1 + (4 — 3/3)(7 — 
+ 12(7 — — 1)], and 7 is the ratio of specific heats (Chen & Taam 1993). 

For the local radiative cooling, we go beyond the model of Abramowicz et al. (1995) 
and adopt the bridging formula introduced by Wandel & Liang (1991, cf. Luo & Liang 
1994). Their formula takes into account free-free and electron scattering opacities, with 
Comptonization, and provides an ad hoc description of the radiative transfer when the 
effective optical depth r e ff ~ 1. A correct description of the r e ff ~ 1 regime requires 
solving the problem of the vertical disk structure. This is notoriously difficult and the 
results depend on the unknown vertical stratification of the viscous heating (e.g. Shaviv 
and Wehrse 1989). Other phenomenological prescriptions for Q- at r e ff ~ 1 (e.g. Lasota 
& Pelat 1991, Narayan & Yi 1995b, Artemova et al. 1994) do not modify the solution 
topology, as we show below. Finally, the equation of state is written as the sum of gas and 
radiation pressure, with p ra( j = Q-(tr/2c) } where tr is the Rosseland mean optical depth. 

While most of the results presented here make use of the above equations, we also 
show for comparison some calculations with a different formalism developed by Narayan 
& Yi (1995b). The major differences in the latter model are (i) the explicit use of a 
self-similar solution with significant sub-Keplerian rotation in the advection-dominated 
regime (Narayan & Yi 1994), (ii) a Newtonian rather than the more accurate pseudo- 
Newtonian potential, (iii) a different vertical integration for the advection term, Q a dv, (iv) 
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the inclusion of magnetic pressure in the equation of state, thermal synchrotron cooling 
and more elaborate Comptonization, (v) allowing the plasma to be two-temperature, and 
(vi) a different form of Q- in the regime r e ff ~ 1. 

Most of these differences lead to only modest changes. However, the inclusion of 
Comptonized synchrotron cooling (item iv), which is a very efficient process at high tem- 
peratures, has a significant effect especially for r e ff < 1. Nevertheless, as discussed below, 
the qualitative form of the results is unaffected compared to the simpler model described 
earlier. The difference (iii) arises primarily because of different assumptions on the profile 
of H(R). The form of £ given in the line below equation (2.4) corresponds to a constant 
H. In the more general case, £ has an additional term —(4 — 3/3)<iln H / din R; setting 
<fln H / din R = 1 recovers the self-similar advection term used in Narayan & Yi (1994). 
Pure gas flows with 7 = 5/3 give non-rotating self-similar behavior according to Narayan 
& Yi (1994), whereas this is not seen in the approach of Abramowicz et al. (1995). 

3. TOPOLOGY OF THERMAL EQUILIBRIA 

The above equations were solved for various choices of the four disk parameters: the mass 
of the central black hole M, the mass accretion rate M, the viscosity coefficient a, and the 
distance from the central object R. It is convenient to discuss the results in the log(M/ME) 
vs. log E plane, where Me = Le/c 2 = (AircGM)/ K es c 2 is the Eddington mass accretion 
rate. (n es is the opacity due to electron scattering.) 

Figure 1 shows sequences of models corresponding to different choices of a. The 
solutions refer to a fixed black hole mass M = 10M©, and to radii R = 5Rq (Fig. la) 
and R = 30Rq (Fig. lb). At each R } there is a particular self-crossing sequence of models 
which corresponds to a critical value a C rit of the viscosity parameter; a C rit is a function of 
R } reaching a minimum value of 0.21 at R ~ 7Rq and reaching unity at R ~ 126 Rq- 

Figure 1 shows that the two a = a C rit curves divide the \og(M / Me) — log S plane 
into four distinct sections, dividing the disk models into four classes. We discuss below 
the properties of these four types of solutions as a function of three parameters, namely, 
the viscosity parameter a (relative to a C rit), the effective optical depth r e ff, and the ratio 
of advective cooling to heating, q = Q&dv/Q+- In the following we use the term slim disk 
to refer to systems with H/R < 1 as well as those with H/R 1. The latter are often 
described as thick disks, but this distinction is not important since the slim disk equations 
appear to describe even thick systems quite well (Narayan & Yi 1995a). 

Type I: Geometrically thin (H/R <C 1) or slim (H/R optically thick (r e fi ^> 1 ) disks 

with small viscosity (a < a C rit )■ The equilibrium sequences are shown by the S curves in 
solid lines at the right of Figs, la, b. Models on the lower and middle branches of the S are 
the usual Shakura-Sunyaev (1973) solutions. Models on the upper branch are dominated 







by radiation pressure, cooled by advection, and have g — > 1 , H/R < 1. This branch was 
discovered by Abramowicz et al. (1988). The lower and upper branches are thermally and 
viscously stable, whereas the middle branch is both thermally and viscously unstable. 

Type II: Geometrically thin (H/R <C 1 ) disks with small advective cooling (q <C 1) and 
large viscosity (a > a C rit )■ These solutions are shown by the dashed lines in Figs, la, b. 
The equilibrium sequence consists of three branches. The two branches on the right are 
the same as the middle and lower branches of Type I (but for a higher value of a) and 
have r e ff > 1. The branch on the left corresponds to a gas pressure dominated, optically 
thin (r e ff < 1), viscously stable, but thermally unstable, flow. Such sequences of models 
were calculated by Shapiro, Lightman & Eardley (1976) and more recently by Luo & Liang 

(1994) . 

Type III: Geometrically slim (H/R< 1 ), optically thin (V e ff <C 1 ) disks with small viscosity 
(a < a C rit )■ These solutions are shown by the dotted lines in Figs, la, b. The equilibrium 
sequence consists of two branches. The branch on the right is similar to the viscously 
stable, but thermally unstable, branch of Type II. The branch on the left is a stable gas 
pressure supported, optically thin, advection-dominated (q — > 1) slim disk (H/R < 1) 
solution. This branch was found independently by Abramowicz et al (1995) and Narayan 
and Yi (1995b) (see also Narayan & Popham 1993). 

Type IV: Geometrically slim (H/R < 1 ) disks with large advection (q — > 1) and large 
viscosity (a > a C rit )■ These solutions, shown by the long-dashed lines in Figs, la, b, are 
fully advection-dominated and extend over the entire range of M / Me , switching from being 
optically thin (r e ff < 1) at low M to optically thick at high M. This sequence corresponds 
to some of the sub-Keplerian, advection-dominated solutions found by Narayan and Yi 
(1994, 1995a), but the fact that they extend without break in M is a new result. 

In order to test whether the results are specific to the model of Abramowicz et al 

(1995) , we have redone the calculations with the model of Narayan & Yi (1995b) briefly 
outlined towards the end of §2. Figure lc shows the topology of the solutions at R = 5Rq 
for a = 0.1 and 1 with this model. We see exactly the same four types of solutions as 
in Figs, la, b. The main differences are that the curves move to lower values of M, 
because of the more efficient cooling due to Comptonized synchrotron emission, while the 
inclusion of magnetic pressure causes minor changes in the geometrically thin disk branch. 
Nevertheless, the qualitative features of the topology are preserved intact, and even the 
value of a cr a is very little changed: 0.41 in Fig. lc against 0.25 in Fig. la. This test shows 
that the classification of solutions described above is quite robust and is independent of 
the details of the particular dynamical or cooling model employed. In fact, even when the 
viscosity coefficient v is taken to be proportional to gas pressure rather than to the total 
pressure, the topology continues to remain the same. All four types of solutions exist, but 
the value of M cr i t increases substantially. 
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4. DISCUSSION OF THE VARIOUS STABLE EQUILIBRIA 

Apart from showing the different topologies, Fig. lc also indicates the stability properties 
of the solutions. The line segments which are shown by dashed lines, and identified by the 
letter U, are thermally unstable and therefore not relevant for real flows. The remaining 
segments, shown by solid lines, are both thermally and viscously stable. Solutions identified 
by the letter C represent cooling-dominated flows. Solutions marked A which occur at 
higher accretion rates describe advection-dominated flows, where cooling is inefficient, 
either because the medium is optically thin or too optically thick. 

For each R and a, the C solutions terminate at a maximum M above which they are 
thermally and viscously (Lightman & Eardley 1974) unstable. The solid lines in Fig. 2 
show this maximum M as a function of R for three choices of a: 0.01, 0.1, 1. There are 
no stable cooling-dominated thin disk solutions above these lines. 

The advection-dominated A solutions come in several varieties. For a < a cr i t} these 
solutions split into two branches, those at small M which are advection-dominated because 
they are optically thin and unable to cool, and those at large M which are extremely 
optically thick and again unable to cool (Narayan & Yi 1994, 1995a). In between these 
two, there is a range of M where there is no advection-dominated solution allowed. The 
new and unexpected feature revealed by the present calculations is that, for a > a cr i t} 
the A solutions span the entire range of accretion rates and are present for all values of 
M. The critical parameter a cr a separates these two topologies and is a function of R. 
The resulting situation is indicated by the dashed and dotted lines in Fig. 2. Advection- 
dominated solutions are allowed only below and to the left of the dashed lines and above 
and to the left of the dotted lines. 

A study of Fig. 2 reveals that four different possibilities arise in the MR plane. We 
identify the corresponding zones as 1-4. In zone 1, the only stable solution allowed is an 
advection-dominated flow, while in zone 2 only the cooling-dominated solution is allowed. 
In zone 3 both solutions are possible. Which of the two does nature pick? Perhaps the 
answer depends on the history of the accreting gas. Alternatively, a thermal instability in 
the upper layers of thin disks may cause such disks to switch to an advection-dominated 
form whenever the latter solution is available (Narayan & Yi 1995b). 

In zone 4, there is no stable solution available. In analogy with models of dwarf- 
nova outbursts (see Cannizzo 1993 for a review) we suggest that in this region of the 
MR diagram, a thermal limit cycle is set up where the flow oscillates between two states 
(Abramowicz et al. 1988). Note, however, that Honma et al. (1991) found that the disk 
does not reach the upper, optically thick branch. 

The exact locations of the various zones in the MR plane depend on the cooling model 
used, as seen from a comparison of Figs, la, c. In the presence of strong Comptonized 
synchrotron cooling, all the lines in Fig. 2 move down by up to 1.5 in log. One result 
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is that systems with M < Me are likely to sample different zones at different radii, so 
that we may expect quite interesting phenomena in these systems accompanied perhaps 
by complex spectral signatures. This may explain the variety of spectral states observed 
in X-ray binaries and active galactic nuclei. 

Figure 2 shows that advection-dominated solutions are available over a surprisingly 
wide range of conditions. Indeed, if the advection-dominated branch is preferred over the 
cooling-dominated one in zone 3, as we speculate above, then it may well be that advection- 
dominated accretion is the most common form of accretion in black hole systems (especially 
at large a), while standard thin disks may occur less frequently. 

We thank Edison Liang for useful discussions and the referee, Paul Wiita, for helpful 
comments. This work was supported in part by NSF Grant AST 9148279. 
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FIGURE CAPTIONS 

Figure 1. — (a) Thermal equilibria of accretion disks at R/Rq = 5. The two heavy 
solid lines correspond to solutions with a = a cr a = 0.2446. The solid, dotted, dashed, 
and long-dashed lines represent respectively the four types of equilibrium disk sequences 
discussed in §3. The calculations were done using the equations given in §2. (b) Same as 
(a) but for R/Rq = 30. (c) Same as (a) but calculated using the model of Narayan & Yi 
(1995b). Note the close similarity to (a), but with all curves shifted down by about 1.5 in 
the log, primarily due to the effect of strong Comptonized synchrotron cooling. 

Figure 2. — Zones in the \og(M / Me) — \og(R/ Rq) plane corresponding to various com- 
binations of stable equilibria, as discussed in §4. 
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